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STOCHASTIC INTEGRATION AND BOUNDARIES 
OF INFINITE NETWORKS 



PALLE E. T. JORGENSEN AND ERIN R J. PEARSE 



Abstract. We study the boundary theory of a connected weighted graph G from the view- 
point of stochastic integration. For the Hilbert space of Dirichlet-finite functions on 
G, we construct a Gel'fand triple S c "Hg Q S'. This yields a probability measure P and 
an isometric embedding of fig, into L-(S',P), and hence gives a concrete representation 
of the boundary as a certain class of "distributions" in S'. In a previous paper, we proved 
a discrete Gauss-Green identity for infinite networks which produces a boundary repre- 
sentation for a harmonic function u: u{x) = EbdG"^i where the sum is understood in 
a liiniting sense. In this paper, we use techniques from stochastic integration to make the 
boundary bdG precise as a measure space, and replace the suin with an integral over S', 
thus obtaining a boundary integral representation for the harmonic function «. 
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13 1. Introduction 

14 In this paper, we develop a boundary theory for an infinite network (connected weighted 

15 graph) G, using some techniques from the theory of stochastic integration. For the Hilbert 

16 space TYg of finite-energy functions on G, we construct a Gel'fand triple So £ 'Wg £ S'^y, 

17 where both containments are strict, and the inclusion mappings are continuous. Here, So 

18 is a space of "test functions" on the network and S'q is a class of "distributions" on the 

19 network, analogous to Schwartz's classical /t/ncf/oni of rapid decay and tempered distri- 

20 butions, respectively. A key result of this paper is Theorem 5.2, which establishes an 
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isometric embedding of "He into the Hilbert space L^{S'q, P), where P is a Gaussian proba- 
biUty measure on S'(^. In Martin boundary theory, elements of the boundary are understood 
as certain minimal harmonic functions. We are studying finite-energy harmonic functions 
instead of positive harmonic functions, but the construction outlined above allows us to 
study elements of the boundary in an analogous fashion. 

In a previous paper, we proved a discrete Gauss-Green identity for infinite networks: 



Here, £ is a Dirichlet form and A is the graph Laplacian; see Theorem 2.21 and the discus- 
sion preceding it for precise definitions of the other terms. Formula (1.1) yields a boundary 
representation for a harmonic function u: 



where C is a constant and the sum is actually defined as a limit of "Riemann sums" over an 
increasing sequence of finite subnetworks of G; see Definition 2.20. In this paper, we use 
functional integration techniques from stochastic integration to make the boundary bd G 
precise as a measure space, and replace the sum with an integral over S'^, thus obtain- 
ing a boundary integral representation (in a sense analogous to that of Poisson or Martin 
boundary theory) for the harmonic function u. 

Another key result of this paper is Corollary 6.1, which follows readily from Theo- 
rem 5.2 and gives a boundary integral representation for harmonic functions of finite en- 
ergy: 



Here, {hx\xeG is a family of harmonic functions parametrized by the vertices and discussed 
in detail just below (and see also Lemma 2.17). Given a transient network, this allows 
one to identify a space of functions in S'q corresponding to the boundary of the network 
(in a manner reminiscent of the Martin boundary). Additionally, Example 7.1 presents the 
construction of a harmonic function of finite energy on a network with one "graph end" 
(in fact, a two-parameter family of such networks). The existence of such functions was 
first proved in [CW92], but we have never seen an exphcit formula given before. We now 
proceed to describe these results in a bit more detail. The reader is also referred to the 
survey paper [JP09e] which gives an overview of how the results of the present paper fit 
into a larger investigation of functions of finite energy on resistance networks, and the 
eff'ective resistance metric. 

1.1. Overview. We define what it means for a function on a network to have finite energy 
in Definition 2.5. Then the (discrete) Dirichlet energy form £ (also given in Definition 2.5) 
is an inner product on the space of functions of finite energy, and in fact produces a Hilbert 
space which we denote by TYe. The space "Kg consists of equivalence classes of functions 
on the vertices of G, where u ^ v iff m - v is a constant function. In a previous paper, 
we constructed a reproducing kernel {vv}j.£co for this Hilbert space and used it to prove a 
discrete Gauss-Green identity which is recalled in Theorem 2.21. 

The space TYfi also enjoys an orthogonal decomposition into the subspace Tin of (£- 
limits of) finitely supported functions and the subspace 'Harm of harmonic functions; see 
Definitions 2.9-2.10 and Theorem 2.12. Since the reproducing kernel behaves well with 




(1.1) 




(1.2) 
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respect to projections, we also have reproducing kernels {/il.iec" for Tin and {/j.ilxeG" for 
Tlarm, where 

fx PnnV.x, and h_, := P-HarmVx- 

It should be noted that these kernels are reproducing up to an additive constant; in other 
words, for some fixed reference vertex o. 



and similarly for /i and /z j. 

Recall the classical result of Poisson that gives a kernel A: : Q x (9Q — » R from which a 
bounded harmonic function can be given via 



We are motivated by the discrete analogue of this result appearing in ( 1 .2). 

Comparison of (1.2) and (1.4) makes one optimistic that bdG can be realized as a 
measure space which supports a measure corresponding to thus replacing the sum in 
(1.2) with a integral. In Corollary 6.1, we do precisely this. 

Boundary theory of harmonic functions can roughly be divided three ways: the bounded 
harmonic functions (Poisson theory), the nonnegative harmonic functions (Martin theory), 
and the finite-energy harmonic functions studied in the present paper. While Poisson theory 
is a subset of Martin theory, the relationship between Martin theory and the study of TYg 
is more subtle. For example, there exist unbounded functions of finite energy; cf. [JP09d, 
Ex. 13.10]. Some results detailing the interrelations are given in [Soa94, §3.7]. 

Whether the focus is on the harmonic functions which are bounded, nonnegative, or 
finite-energy, the goals of the associated boundary theory are essentially the same: 

(1) Construct a space D which extends the original domain T)\ this can be done by 
taking closure, compactifying, or similar operations. 

(2) One can then identify the boundary bd £) as £) \ 2), or (if the boundary thus ob- 
tained would be larger than necessary/practical for the application in mind), as 
some subset of T)\D. 

(3) Define a procedure for extending harmonic functions u from £) to bdD. This 
extension u may be a measure (or other linear functional) on bd £); it may not be 
a function. 

(4) Obtain a kernel ]k;(x,y6) defined on D x bd£) against which one can integrate the 
extension u so as to recover the value of m at a point in 2): 



whenever m is a harmonic functions of the given class. 

Our approach to (1) is to use Gel'fand triples to extend the original domain, a method 
which is novel as far as we know. In a forthcoming work [JPIO], we will introduce an 
interpolation formula that uses the analytic framework developed in this paper, and which 
turns G into a stochastic process. In particular, the interpolation formula allows one to find 
continua which naturally extend G. 



{Vx, w)£ = U{x) - u(o), 




(1.4) 
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The difference between our boundary theory and that of Poisson and Martin is rooted 
in our focus on "He rather than both classical theories concern harmonic functions with 
growth/decay restrictions. By contrast, provided they neither grow too wildly nor oscillate 
too wildly, elements of may have values tending to both +00 and -00. See [JP09d, 
Ex. 13.10] for a function h e 'Harm which is unbounded in this way. Positive harmonic 
functions are naturally given to analysis based on probabilistic and potential-theoretic 
techniques, and the companion study of superharmonic (or subharmonic) functions is in- 
dispensable. Without positivity, however, one can get more mileage by considering the 
Dirichlet form & as an inner product and studying the resulting Hilbert space geometry. 

Our boundary essentially consists of (equivalence classes of) infinite paths which can 
be distinguished by monopoles, i.e., two paths are not equivalent iff there is a monopole 
w with different limiting values along each path. It is an immediate consequence that 
recurrent networks have no boundary, and transient networks with no nontrivial harmonic 
functions have exactly one boundary point (corresponding to the fact that the monopole 
at X is unique). In particular, the integer lattices (Z'', 1) each have 1 boundary point for 
d > 3 and boundary points for d - 1,2. In contrast, the Martin boundary of (Z'', 1) is 
homeomorphic to the unit sphere 5" ' (where 5" = {-1, 1)), and each (Z'', 1) has only one 
graph ends (except for (Z, 1), which has two); cf. [PW90, §3.B], for example. 

1.2. Outline. In our version of the program outlined above, we follow the steps in the or- 
der (2)-(3)-(l). A brief summary is given here; further introductory material and technical 
details appear at the beginning of each subsection. 

§2 recalls basic definitions and some previously obtained results. 

§3 describes two methods for constructing a Gel'fand triple. The technique presented 
in §3. 1 works for any network (G, c) and makes use of an orthonormal basis of "Kg derived 
from the energy kernel {vil.vec via the Gram-Schmidt algorithm, or equivalently, from the 
domain of a certain operator N. The approach given in §3.2 works only for networks 
where A is an unbounded operator on "Kg. This version of So is constructed in terms of 
the domain of A. 

§4 studies the structure of Sa (the space of test functions) and S'q (the space of distri- 
butions) and establishes some key lemmas for later use. 

§5 proves a key result: Theorem 5.2, which establishes the isometric embedding of 'Hs 
into L~(S'^, P) given by the Wiener transform. Applying this isometry to the energy kernel 
{vj}, we get a reproducing kernel ]k;(x, d¥) given in terms of a version of Wiener measure. 
In fact, P is a Gaussian probability measure on S'q whose support is disjoint from Tin. 
The results in this section hold for any Gel'fand triple; in particular, for either of the ones 
constructed in §3 and §3.2. 

§6 We consider certain measures fix, defined in terms of the kernel and the Wiener 
measure just introduced, which are supported on S^ITin and indexed by the vertices x e 
G". Then points of bdG correspond to limits of sequences [fix,,] where x„ 00, modulo a 
suitable equivalence relation. 

Remark 1.1. While Doob's martingale theory works well for harmonic functions in L°° or 
L^, the situation for TYg is different. The primary reason is that 'Kg is not immediatelly real- 
izable as an space. A considerable advantage of our Gel'fand-Wiener-Ito construction is 
that TYg is isometrically embedded into L^{S'q, P) in a particularly nice way: it corresponds 
to the polynomials of degree 1. See Remark 5.5. 

Another contrast is that A may, in general, be unbounded in our context. Recall that 
when studying an operator, an important subtlety is that "the" adjoint A* depends on the 
choice of domain, i.e., the linear subspace dom(A) c 'H. 
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1 Boundary theory is a well-established subject; the deep connections between harmonic 

2 analysis, probability, and potential theory have led to several notions of boundary and we 

3 will not attempt to give complete references. However, we recommend [Saw97] for intro- 

4 ductory material on Martin boundary and [DS84, LP09] for introductory material on resis- 

5 tance networks. Additionally, [Lyo83,Car73, WoeOO], and the foundational paper [NW59] 

6 provide more specific background. With regard to infinite graphs and finite-energy func- 

7 tions, see [Soa94,SW91,CW92,Dod06,PW90,PW88, Woe86,Tho90]. But we ask different 

8 questions here, and the operator theory we use is different; it does not easily compare to 

9 earlier literature. For some recent related areas, see e.g., [AL08, AAL08, AD06] repro- 

10 ducing kernels, [Arv86] Markov operators, [Cho08] graph analysis, and [AP09] operator 

1 1 theory. 

12 There has been a recent interest in analysis and potential theory on infinite-dimensional 

13 spaces, and the use of stochastic integration in conjunction with reproducing kernels [HNS09, 

14 XialO, CdVTUlO, ZXZ09], and Gel'fand triples [HLWIO, BKO07]. Although our setting 

15 here is different, we are able to adapt these tools for the task at hand. This is nontriv- 

16 ial because, in the classical case, there is a natural differentiable structure around, and 

17 therefore the choice of a Schwartz space going into a useful Gel'fand triple is often rather 

18 conventional. But by contrast, we deal with discrete structures, and so we must give up 

19 differential operators. Nonetheless, we exhibit Schwartz spaces that yield Gel'fand triples 

20 which accomplish what we need. 



21 2. Basic terms and previous results 

22 We now proceed to introduce the key notions used throughout this paper: resistance 

23 networks, the energy form fi, the Laplace operator A, and their elementary properties. Our 

24 approach is somewhat different from existing studies of networks in the literature, and 

25 so we take this opportunity to introduce the tools we will need: an unbounded Laplace 

26 operator with dense domain in a Hilbert space, a two-point reproducing kernel for this 

27 Hilbert space, the quadratic form associated to the Laplacian, and Gelfand triples. Since 

28 these are tools not commonly used in geometric analysis, we include their definitions and 

29 some theorems from earlier papers which will be needed later Additionally, we will use 

30 the theorems of Bochner (Theorem 2.23), and Minlos (Theorem 2.24). 

31 Definition 2.1. A resistance network is a connected graph [G, c), where G is a graph with 

32 vertex set G°, and c is the conductance function which defines adjacency by x ~ y iff 

33 c„ > 0, for x,y e G". We assume c^ - Cy.v 6 [0, oo), and write c{x) := Yjv-x'^-^y 

34 We require c{x) < oo (note that we allow vertices of infinite degree), but c{x) need not 

35 be a bounded function on G". The notation c may be used to indicate the multiplication 

36 operator (cv)(x) := c{x)\'{x), i.e., the diagonal matrix with entries c{x) with respect to the 

37 (vector space) basis 

38 In this definition, connected means simply that for any x,y E G*^, there is a finite se- 

39 quence {j^/I'Lq with x - xq, y - x„, and Ct,_,.i, > 0, / = !,...,«. Conductance is the 

40 reciprocal of resistance, so one can think of (G, c) as a network of nodes connected by 

41 resistors of resistance c^^ . We may assume there is at most one edge from x to y, as two 

42 conductors c\ and ci, connected in parallel can be replaced by a single conductor with 



6 



PALLE E. T. JORGENSEN AND ERIN R J. PEARSE 



conductance c^y - c',. + Cj.,,. Also, we assume c.o = so that no vertex has a loop, as 
electric current will never flow along a conductor connecting a node to itself. ' 

Definition 2.2. The Laplacian on G is the linear diff'erence operator which acts on a func- 
tion V : G" R by 

(Av)(x):=^c,,(v(x)-v(3;)). (2.1) 

A function V : G° R is harmonic iff Av(x) = for each x € G . 

We have adopted the physics convention (so that the spectrum is nonnegative) and thus 
our Laplacian is the negative of the one commonly found in the PDE literature. The net- 
work Laplacian (2. 1) should not be confused with the stochastically renormalized Laplace 
operator c 'A which appears in the probabiUty literature, or with the spectrally renormal- 
ized Laplace operator c '^^Ac '^^ which appears in the literature on spectral graph theory 
(e.g., [ChuOl]). 

Definition 2.3. An exhaustion of G is an increasing sequence of finite and connected sub- 
graphs {Gtl^j, so that Gk c G^+i and G = IJGk- Since any vertex or edge is eventually 
contained in some Gk, there is no loss of generality in assuming they are contained in Gi, 
for the purposes of a specific computation. 

Definition 2.4. The notation 

A-eG» -teG,- 

is used whenever the limit is independent of the choice of exhaustion {Gk} of G. This is 
clearly justified, for example, whenever the sum has only finitely many nonzero terms, or 
is absolutely convergent as in the definition of S in Definition 2.5. 

Definition 2.5. The energy of functions m, v : G" ^ C is given by the (closed, bilinear) 
Dirichlet form 

6(u, V) 5 2 Z ^■^■J'^"^^^ " ^(y))(v(x) - v(y)), (2.3) 

xeC° veG" 

with the energy of u given by £(«) :- &{u, u). The domain of the energy is 

dom£={M:G°^C;£(M)<oo). (2.4) 

Since c„ - Cy^ and c„ = for nonadjacent vertices, the initial factor of i in (2.3) 
implies there is exactly one term in the sum for each edge in the network. 

2. 1 . Tlie energy space "Kg. Let 1 denote the constant function with value 1 and recall that 
ker£ = CI. 



Nonetheless, self-loops may be useful for technical considerations: one can remove the periodicity of a 
random walk by allowing self-loops. This can allow one to obtain a "lazy walk" which is ergodic, and hence 
more tractable. See, for example, [LPW08,LP09]. 
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Definition 2.6. The energy form £ is symmetric and positive definite on dom£. Then 
dom£/Cl is a vector space with inner product and corresponding norm given by 

<M, v}e := £(m, v) and ||m||£ := £(m, u)^'^. (2.5) 
The energy Hilbert space "He is dom £/Cl with inner product (2.5). 

Definition 2.7. Let v i be defined to be the unique element of fie for which 

{vx, m)£ = u{x) - u{o), for every u e 'He. (2.6) 
The collection {Vx}jcec'' forms a reproducing kernel for 'Kg ( [JP09b, Cor. 2.7]); we call it 
the energy kernel and (2.6) shows its span is dense in "Hs- Note that Vo corresponds to a 
constant function, since {vo, u}s - for every u e 'He. Therefore, v„ is often ignored or 
omitted. 

Definition 2.8. A dipole is any v e •Kg satisfying the pointwise identity Av - 5x - Sy for 

some vertices x,y e One can check that Avx = 6x - 6,,', cf. [JP09b, Lemma 2.13]. 

Definition 2.9. For v e 'He, one says that v has finite support iff there is a finite set F QG^ 
for which v{x) - k eC for all x i F, i.e., the set of functions of finite support in "He is 

span{(5^) = {m e dom£ : u{x) - k for some k, for all but finitely many x e G"}, (2.7) 

where Sx is the Dirac mass at x, i.e., the element of "Kg containing the characteristic func- 
tion of the singleton [x]. It is immediate from (2.3) that &{5x) - c{x), whence 6x e "He. 
Define Tin to be the closure of span{(5^ ) with respect to £. 

Definition 2.10. The set of harmonic functions of finite energy is denoted 

-Harm {v € 'He '■ Av(jc) = 0, for all x e G"). (2.8) 
Note that this is independent of choice of representative for v in virtue of (2.1). 

Lemma 2.11 ( [JP09b, 2.11]). For any x e G", one has {6x, u}e = Am(-«). 

The following result follows easily from Lemma 2. 1 1; cf. [JP09b, Thm. 2. 15]. 

Theorem 2.12 (Royden decomposition). 'He - Tin © 'Harm. 

Definition 2.13. A monopole at jc e G° is an element Wi e 'He which satisfies Awxiy) - 
6xy, where 6xy is Kronecker's delta. In case the network supports monopoles, let Wo always 
denote the unique energy-minimizing monopole at the origin. 

When "Hg contains monopoles, let Mx denote the vector space spanned by the monopoles 
at X. This implies that Mx may contain harmonic functions; see [JP09b, Lemma 4.1]. With 
Vi and fx = PrinVx as above, we indicate the distinguished monopoles 

w'^, := v.v + Wo and Wx :- fx + Wo. (2.9) 

Remark 2.14. Note that w,, e Tm, whenever it is present in 'He, and similarly that w'x is the 
energy-minimizing element of Mx. To see this, suppose Wx is any monopole at x. Since 
Wx e 'He, write Wx = f + lihy Theorem 2.12, and get S(wx) - 8(f) + &{h). Projecting 
away the harmonic component will not affect the monopole property, so w^x - PrmWx is 
the unique monopole of minimal energy. The Green function is g{x,y) - w°{x), where w° 
is the representative of Wj, satisfying w°(o) - 0. 
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Definition 2.15. The dense subspace of He, spanned by monopoles (or dipoles) is 

M := span{v J ,£G0 + span{w^, w^tlAec" ■ (2.10) 

Let Ayvi be the closure of the Laplacian when taken to have the dense domain M. Since 
A agrees with A^^ pointwise, we may suppress reference to the domain for ease of notation. 

Lemma 2.16 ( [JP09b, Lemma 3.5]). Am is Hermitian with (u, Amu)s > 0/or all u e Ai. 

Lemma 2.17 ( [JP09b, Lemma 3.6]). When the network is transient, M contains the 
spaces span{vv), span{/t}, anc/ span{/zr), where = PrmVx and hx — P-HannVx- When the 
network is not transient, M — span{v.v) = spanf/j^}. 

Remark 2.18 (Monopoles and transience). The presence of monopoles in 'He is equivalent 
to the transience of the simple random walk on the network with transition probabilities 
p{x,y) = Cxy/cix): note that if Wx is a monopole, then the current induced by Wx is a unit 
flow to infinity with finite energy. It was proved in [Lyo83] that the network is transient if 
and only if there exists a unit current flow to infinity; see also [LP09, Thm. 2.10]. 

2.2. The discrete Gauss-Green identity. The space M is introduced as a dense domain 
for A and as the scope of validity for the discrete Gauss-Green identity of Theorem 2.21. 

Definition 2.19. If // is a subgraph of G, then the boundary of H is 

hdH -.^ {x e H - ly e H^,y ~ x). (2.11) 
The interior of a subgraph H consists of the vertices in H whose neighbours also lie in H: 

intH -.^ {x e H y ~ X => y e H} ^ H\hdH. (2.12) 
For vertices in the boundary of a subgraph, the normal derivative of v is 

l^(x) := Yj '^-^(^W - ^W), for xehdH. (2. 13) 

yeH 

Thus, the normal derivative of v is computed like Av(ji:), except that the sum extends only 
over the neighbours of x which lie in H. 

Definition 2.19 will be used primarily for subgraphs that form an exhaustion of G, in 
the sense of Definition 2.3: an increasing sequence of finite and connected subgraphs [Gk], 
so that Gk c Gk+i and G = \JGif Also, recall that 2bdc lini/t^oo SbdCi from Defini- 
tion 2.20. 

Definition 2.20. A boundary sum is computed in terms of an exhaustion {Gk] by 

bdC bdCj 

whenever the limit is independent of the choice of exhaustion, as in Definition 2.4. 

On a finite network, all harmonic functions of finite energy are constant, so that 'He - 
Tin by Theorem 2.12, and one has £(m, v) = Yjxec" u{x)Av(x), for all u,v e He- In fact, 
this remains true for recurrent infinite networks, as shown in [JP09b, Thm. 4.4]; see also 
[KY89]. However, the possibilities are much richer on an infinite network, as evinced by 
the following theorem. 
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Theorem 2.21 (Discrete Gauss-Green identity). Ifue Us and v e At, then 

(m,v)£ = 2mAv + 2m£. (2.15) 

C" bdC 

The discrete Gauss-Green formula (2.15) is the main result of [JP09b]; that paper con- 
tains several consequences of the formula, especially as pertains to transience. 

2.3. Gel'fand triples and duality. One would like to obtain a (probability) measure space 
to serve as the boundary of G. It is shown in [Gro67, Gro70, Min63] that no Hilbert space 
of functions is sufficient to support a Gaussian measure P (i.e., it is not possible to have 
< PCH) < oo for a cr-finite measure). However, it is possible to construct a Gel'fand 
triple (also called a rigged Hilbert space): a dense subspace SofH with 

SQ'HqS', (2.16) 

where S is dense in "K and S ' is the dual of S . Additionally, S and S ' must also satisfy 
some technical conditions: 5 is a dense subspace of with respect to the Hilbert norm, 
but also comes equipped with a strictly finer "test function" topology. When 5 is a Frechet 
space equipped with a countable system of seminorms (stronger than the norm on Iri), then 
the inclusion map of S into is continuous; in fact, it is possible to chose the seminorms 
in such a way that one gets a nuclear embedding (details below). Therefore, when the dual 
S ' is taken with respect to this finer (Frechet) topology, one obtains a strict containment 
'HQS'. It turns out that 5" is large enough to support a nice probability measure, even 
though 7f is not. 

It was Gel'fand's idea to formalize this construction abstractly using a system of nu- 
clearity axioms [GMS58,Min58,Min59]. Our presentation here is adapted from quantum 
mechanics and the goal is to realize bdG as a subset of 5". We will give a "test function 
topology" as a Frechet topology defined via a specific sequence of seminorms, using either 
an onb for 'He (in §3) or the domain of A'^ (in §3.2). 

Remark 2.22 (Tempered distributions and the Laplacian). There is a concrete situation 
when the Gel'fand triple construction is especially natural: 'H - L^{R,dx) and S is the 
Schwartz space of functions of rapid decay. That is, each / e 5 is C°° smooth functions 
which decays (along with all its derivatives) faster than any polynomial. In this case, S is 
the space of tempered distributions and the seminorms defining the Frechet topology on S 
are 

p,„(f) := sup{|//"Hx)| ■■.xeR,0<k,n< m), m = 0,1,2,..., 

where /'"^ is the n* derivative of /. Then S ' is the dual of S with respect to this Frechet 
topology. One can equivalently express S as 

5 := {/ 6 \ (p2 + e^)"/ e V«), (2.17) 

where P : f{x) ^ and Q : f(x) xf(x) are Heisenberg's operators. The operator 
+ Q- is often called the quantum mechanical Hamiltonian, but some others (e.g., Hida, 
Gross) would call it a Laplacian, and this perspective tightens the analogy with the present 
study. In this sense, (2.17) could be rewritten S :- domA°°; compare to (3.15). We 
show that a general network (G, c) always has a harmonic oscillator; in fact, we discuss 
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an operator N in Definition 3.7 which is unitarily equivalent to P- + and hence has the 
same spectrum. 

The duality between S and S ' allows for the extension of the inner product on "TY to a 
pairing of S and S': 

(-,■>/,: X -H ^ C to {;■}:„: S xS' ^R. 

In other words, one obtains a Fourier-type duality restricted to S . Moreover, it is possible 
to construct a Gel'fand triple in such a way that P(5'') = 1 for a Gaussian probability 
measure P. When applied to •?/ = "TYg, the construction yields three main outcomes: 

(1) The next best thing to a Fourier transform for an arbitrary graph. 

(2) A concrete representation of 'T/g as an measure space "TYg = L^iS^F). 

(3) A boundary integral representation for harmonic functions of finite energy. 

As a prelude, we begin with Bochner's Theorem, which characterizes the Fourier trans- 
form of a positive finite Borel measure on the real line. The reader may find [RS75] helpful 
for further information. 

Theorem 2.23 (Bochner). Let G be a locally compact abelian group. Then there is a 
bijective correspondence T '■ M(G) — > PD{G), where Ai(G) is the collection of measures 
on G, and PD{G) is the set of positive definite functions on the dual group ofG. Moreover, 
this bijection is given by the Fourier transform 

r-.v^ify by ifyi^)^ f e'<^-'^ dv{x). (2.18) 

Jc 

For our representation of the energy Hilbert space in the case of general resistance 
network, we will need Minlos' generalization of Bochner's theorem from [Min63, Sch73]. 
This important result states that a cylindrical measure on the dual of a nuclear space is 
a Radon measure iff its Fourier transform is continuous. In this context, however, the 
notion of Fourier transform is infinite-dimensional, and is dealt with by the introduction of 
Gel'fand triples [Lee96]. 

Theorem 2.24 (Minlos). Given a Gelfand triple S Q 'H Q S', Bochner's Theorem may 
be extended to yield a bijective correspondence between the positive definite functions on 
S and the Radon probability measures on S '. Moreover, in a specific case, this correspon- 
dence is uniquely determined by the identity 

ei<".f>'H dF(^) = e-iM'H^ (2.19) 

where (■, is the original inner product on 'H and (■, ■),^ is its extension to the pairing 
onS xS'. 

Formula (2.19) may be interpreted as defining the Fourier transform of P; the function 
on the right-hand side is positive definite and plays a special role in stochastic integration, 
and its use in quantization. 

3. Gel'fand TRIPLES FOR 'Kg 

In this section, we describe two methods for construction a Gel'fand triple for 'He- The 
first method is applicable to all networks, but relies on the choice of some enumeration 
of the vertices of G, and the Gram-Schmidt algorithm for producing an onb. However, 
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we will see that the Gram-Schmidt algorithm yields a much more explicit formula than 
usual, in the present context. The second method is applicable only when the Laplacian is 
unbounded. However, in this case the construction does not require any enumeration (or 
onb) and may provide for more feasible computations. 

Remark 3.1. Note that Sa and 5^ consist of R- valued functions in this section. This tech- 
nical detail is important because we do not expect the integral J^, from (2.19) to 
converge unless it is certain that {u, ■) is R-valued. After the Wiener embedding is carried 
out in Theorem 5.2, all results can be complexified. 

3.1. Gel'fand triples via Gram-Schmidt. In this section, we describe a Gel'fand triple 
for "He where the class of test functions Sq is described in terms of the decay properties 
of a certain orthonormal basis (onb) for TYg. We will see in Remark 5.9 that this onb 
corresponds to a system of i.i.d. random variables (which are, in fact, Gaussian with mean 
and variance 1). 

The onb {e„)„eN comes by applying the Gram-Schmidt process to the reproducing kernel 
{i'A„);ieN, where we have fixed some enumeration {x„)„eH of the vertices G \ {o}. That is, we 
put xq - o and henceforth exclude xq from the discussion, as it will not be relevant. Given 
{ei , . . . , e„-i ), one obtains e„ via 







\ Me' 








. 







V,-, 






M2,i 


M2,2 













V'.v, 


£3 




M3.1 


M3,2 


M3_3 


















M„,2 


M„,3 








- ^'-r,, - 



(3.1) 



Consequently, for each e N, the triangular nature of M gives 



spanK,, . . .,VxJ = span{ei, . . .,5^}. (3.2) 

Remark 3.2. Note that the reproducing kernel gives one an explicit formula for the entries 
of the inverse of this particular Gram-Schmidt matrix: 

(M-% = (v,„ = sjixd - sj(o). (3.3) 

This is certainly in distinct contrast with the general case, and allows us to find a formula 
for the entries of M itself in Lemma 3.3. 

Lemma 3.3. The entries of the Gram-Schmidt matrix M are given by 

{{Aedixj), j <i 
Mij^\ for I, J ^ 1,2,.... (3.4) 

10 else. 

Proof. For j < i, an application of (2.11) gives 

Asiixj) = (5,,, e,)e = U.v,, ^ M,-,iV.,J = ^ M,;t(5,. , v.,,)6 = ^ Ma(5x;ta) - ^x/o)), 

^ k<i ' g k<i k<i 

where the last equality comes by (2.6). Note that d^iy) = 6x{o) for every y except y - x,so 
the last sum above has a nonzero term only for k = j, and the result follows. □ 
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From (3.3) and Lemma 3.3, we have the handy conversion formulas: 



Si - ^ Ae,(xj)v.v. and v,,. - ^ {ek(_Xi) - etio)) sj 



Lemma 3.4. We have the identity 



^ (skixi) - Skio)) Aek(xj) = 6ij, for ; = 1, 2, . . . . 

j<k<i 

Proof. By formula (3.3), the left side of (3.6) is equal to 



(3.5) 



(3.6) 



^ isk(xi) - Ek{o))Ask(xj) = A '^{v_^.,Sj}Eek(xj) 

j<k<i 



\ k<i 



= Av.iX-Xj) = (5a,(^>) - 5o(_Xj). 



Note that Askixj) - for j > k, so the second sum runs over all k < i. Also, note that 
Sxii-^j) " So{xj) - 6ij for /, y > (and the indexing of M begins at 1, not 0). □ 

Lemma 3.4 can also be proven by combining the identities in (3.5). 

Lemma 3.5. Let Vv.v := {vx,Vy}s, and let E - M"' be defined as in (3.3). Then EE* = V. 

Proof. Computing entrywise, 

(EE*)ij = ^ E^.^^^E_,.^^^ = ^ [ek(xi) - Ek{o)) (et(xj) - Sk(o)) = ^(v^:,, ei>£<v.t,, Sk}E, 

k k k 

which is equal to (vv,, v.v.)fi by Parseval's identity. □ 

Definition 3.6. The space of test functions and the space of distributions corresponding to 
the onb {e„}„eN are defined by 



5g = {i = ^ i„e„ ; V/7 6 N, 3C > such that |i„| < C/«''}, and (3.7) 

neN 

S'g^{^^Yj ^"«« : 3p e N, 3C > such that \^„\ < CnP}. (3.8) 

neN 

Thus, So - {^pen{^ '■ \\A\p < °°] where the Frechet p-seminorm of i = 2„eN ■^nSn is 

n1/2 



Y,'l"\Sn\ 



s € Sg,p e N. 



(3.9) 



Note that the system of seminorms (3.9) is equivalent to system of seminorms defined by 



\\s\\p sup«' \Sn\.. 

neN 



s eSG,p e . 



(3.10) 



in the sense that both define the same Frechet topology on Sa- (Each seminorm in one 
system is dominated by one from the other, but with a different p.) We occasionally find it 
more convenient to calculate with (3.10) instead of (3.9). 
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1 Definition 3.7. Let :- spanlvil vec and define a mapping A/' : — > "Kg by 



n 



kekix„)ek. 



(3.11) 



2 Remark 3.8. From (3.11), one has 



\\Nv,J\l^Y,'^\s,(x„)\' and 



k=l 



n 




3 Note that Sk € by (3.2), and that Nek - ks^ for each ^ € N. We use the symbol N for 

4 the operator discussed in this section by way of analogy with the number operator N from 

5 quantum mechanics. Indeed, N can also be defined as a*a for a suitable operator a and its 

6 adjoint. 

7 In the following lemma, we use the symbol N to denote the closure of the operator N 

8 (i.e., the domain is the closure of span{£„} with respect to the graph norm). 

9 Lemma 3.9. The mapping N is essentially self-adjoint, and is unbounded if and only ifG 

10 is infinite. Moreover, if we define the seminorms p„iu) :— \\iN)"u\\E, then {p„] and {|| • \\p] 

11 induce equivalent topologies on Sg, so that 



12 and u e Sa if and only if p„{u) < oo for each n 6 N. 

13 Proof. Unitary equivalence of "TYg with ^^(Z+) is given by [/ : e„ i-> S,„ where 5„{m) :- 

14 S„j„ (Kronecker <5) for Define : fi(Z+) fi(Z+) by N+6„ = nd„ so that 

15 N+U = [/A/' holds on the dense subspace span{e„). The rest follows by [Sim79, §2]. □ 

16 Corollary 3.10. The inclusion mapping So ^ 'Wg is nuclear, and so Sg £ 'He c S'q is a 

17 Gel fand triple. 

IB Proof. When the space of test functions is defined as domr°° for some operator T with 

19 pure point spectrum (as in (3. 13)), then nuclearity follows if there is a /9 e Z+ such that the 

20 reciprocal eigenvalues of T are p-summable; see [Sim79]. Since N has spectrum Z+ and 

21 Yj7=i " < °° f*^'" P ^ 2, the conclusion follows. □ 

22 Lemma 3.11. The energy kernel {Vx]xeco Frechet-dense subset of Sq. 

23 Proof. In the expansion with respect to the onb as in (3.7), the basis element st has co- 

24 efficients i„ = dnt- Since this sequence {.?„) vanishes after « - k, it clearly satisfies the 

25 required decay condition \sn\ < Cn^'' . From (3.2), the same clearly holds for Vx^. This 

26 shows that the kernel is contained in Sq- 

27 To see that {vv} is dense in Sa, it suffices by (3.2) to show that the onb {e„)„eN is dense. 

28 Given any m = 2 6 Sa and any p e N, there is a C = Cp such that \uk\ < Clk''*^ . Now 

29 if - Y!k=i ^'kSk is the A^* truncation of u, then 




(3.13) 



II" - /lip = supF|M;t -/*| = supFlMt -/tl < kP—j 



k k>N 



^ 0. 



30 Thus, one can always approximate u e Sq by Uf^ 6 span{ei, . . . , e^i). 



□ 
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3.2. Gel'fand triples in the case when A is unbounded. In the case when A : "TYg — » 
'He is unbounded, there is an alternative construction of Sa and S'^^, which begins by 
identifying a certain subspace of A1 = dom Am (as given in Definition 2.15) to act as the 
space of test functions. 

Definition 3.12. Let A-y be a self-adjoint extension of Ayvi; since Am is Hermitian and 
commutes with conjugation (since c is R-valued), a theorem of von Neumann's states that 
such an extension exists. 

Let A^M := (A'vA'v . . . A'v)« be the p-fold product of A^' applied to u e "Kg. Define 
dom(A^) inductively by 

dom(A^) {u ■. 6 dom(A<v)). (3.14) 

Definition 3.13. The (Schwartz) space of potentials of rapid decay is 

Sa-.^domi^), (3.15) 

where dom(A^) :- 0)^=1 dom(AfJ,) consists of all R-valued functions u e "TYg for which 
A^M 6 "Hg for any p. The space of Schwartz distributions or tempered distributions is the 
dual space <Sq of R-valued continuous linear functionals on So- 

Remark 3.14. Note that So is dense in dom(A'v) with respect to the graph norm, by standard 
spectral theory. For each p e N, there is a seminorm on So defined by 

IImII^ I^Mllg. (3.16) 

Since (dom A^, || • ||p) is a Hilbert space for each p eN, the subspace So is a Frechet space. 
Note also that since A is unbounded. So is a proper subspace of "TYg. 

Lemma 3.15. If deg{x) is finite for each x € G^, then one has Vj- e So- 
Proof If deg(x) < oo then Lemma 3.16 shows that 6x e span{vv}. □ 

Note that we take deg(x) < oo as a blanket assumption, in fact, as part of the definition 
of a network. However, it is emphasized in Lemma 3.15 because this is the only place 
where it is really necessary. (But note that deg(x) may be unbounded.) 

Lemma 3.16. For any x e G'^, 6x - c(x)vx - Yjy~x<^xvVy- 

Proof Lemma 2.1 1 implies {5x, m)£ = {c(x)vx - Yj\~x ^xvVy, «)£ for every u e "Hg, so apply 
this to M = Vj, z € G°. Since 6x, Vx 6 "Hg, it must also be that Yjy~x CxyVy e "TYg. □ 

Remark 3.17. When the hypotheses of Lemma 3.15 are satisfied, it should be noted that 
span{v.v) is dense in So with respect to fi, but not with respect to the Frechet topology in- 
duced by the seminorms (3.16), nor with respect to the graph norm. One has the inclusions 



Vx 

Aa(V'.i 



s 

A-vs 



(3.17) 

where s e So and u E TYg. The second inclusion is dense but the first is not. 



u 

A^'M 
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4. The structure of Sa and S'q 

From this point on, we assume that a Gel'fand triple has been chosen, using either 
of the methods described in the previous section. Henceforth, we use the symbol A to 
denote the operator N = N* or the operator At/, depending on how the Gel'fand triple was 
constructed: 



[N, Definitions.? 
A:=i (4.1) 
[At., Definition 3.12. 

4. 1 . The structure of So- We establish that So is a dense analytic subset of "He, and that 
the energy product can be extends not just to a pairing on So x S'q, but all the way to a 
pairing on "TYg x S'^j. Parts of this subsection closely parallel the general theory, and a good 
reference would be [Sim79] or [HidSO]. 

Definition 4.1. Let ^'[a h] denote the usual indicator function of the interval [a, b] c R, and 
let S be the spectral transform in the spectral representation of A, and let E be the asso- 
ciated projection- valued measure. Then define E„ to be the spectral truncation operator 
acting on Tfe by 



Jlln 



E„u := Q*Xfi „iSm - \ E{dt)u. 

I In 

Lemma 4.2. With respect to &, Sa is a dense analytic subspace ofHs- 

Proof. This essentially follows immediately once it is clear that E,, maps "TYg into So- For 
M € 'Hs, and for any p = 1,2,..., 

\\APE,M\l^ r A^"\\E(dA)u\\l < n^"\\u\\l, (4.2) 

Jl/n 

So E„u E Sa- It follows that \\u - E„u\\g, ^ by standard spectral theory. □ 

Tlieorem 4.3. Sa £ 'J^e C >Sq is a Gel'fand triple, and the energy form {■, ■)£ extends to a 
pairing on Sa X >Sq defined by 

{u,vU,:^{KPu,K-Pv)e, (4.3) 
where p is any integer such that \v{u)\ < A'||A''m||£ /or all u E Sa, for some K > 0. 

Proof. In combination with (3.15)-(3.16), Lemma 4.2 establishes that Sa C 'T/g Q S'q is 
a Gel'fand triple. If v £ S'q, then there is a C and p such that Ks, v)tv| < CHA'^sUfi for 
all s E Sq. Set (p{A''s) := {s,v}n, to obtain a continuous linear functional on "Hg (after 
extending to the orthogonal complement of span{A''i} by if necessary). Now Riesz's 
lemma gives a w E "Kg for which {s,v}'w = {APs,w)e for all s e Sa and we define 
A^Pv := w E TYg to make the meaning of the right-hand side of (4.3) clear. □ 

Lemma 4.4. The pairing on Sa x S'q is equivalently given by 

<M,^W= limf(£„M), (4.4) 
where the limit is taken in the topology of S'^^. Moreover, u(^) — {u,^)^ is 'K-valued on S^. 
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Proof. En commutes with A. This is a standard result in spectral theory, as E„ and A 
are unitarily equivalent to the two commuting operations of truncation and multiplication, 
respectively. Therefore, 

Standard spectral theory also gives EnV — » v in 'He,, so 

lim ^{E„u) = \m\{KPs,E„K-''^)E = <A''m, A"''v)£. 

Note that the pairing (■ , is a limit of real numbers, and hence is real. □ 

Corollary 4.5. E„ extends to a mapping E„ : S'q — > 'He defined via {ii,E„^)s :— ^(E„u). 
Thus, we have a pointwise extension of {■ , ■)'w to "TYg X S'q given by 

lim(M,£„^)£. (4.5) 

4.2. The structure of S'q. The next results are structure theorems akin to those found 
in the classical theory of distributions; see [Str03, §6.3] or [AG92, §3.5]. If He, c Sq, 
then Theorem 4.6 would say S'q - \Jp^''{'He) (of course, this is typically false when 
'Harm it 0). 

Theorem 4.6. The distribution space S'q is 



S'q = {^(m) = {A''u, v)e ■■ 3v e'He,pe Z+, Vm e Sq]. (4.6) 

Proof. It is clear from the Schwarz inequality that ^(«) = {A'^u, v)e defines a continuous 
linear functional on Sq, for any v E 'He and nonnegative integer p. For the other direction, 
we use the same technique as in Lemma 4.3. Observe that if ^ E S'q, then there exists K, p 
such that |^(m)| < /ir||A''M||£ for every u E Sq. This implies that the map ^ : K^u i-> ^(m) is 
continuous on the subspace Y - span{A^M : u E 'He, p E Z^). This can be extended to all 
of 'He by precomposing with the orthogonal projection to Y. Now Riesz's lemma gives a 
V E 'He for which ^(m) = {K^u, v)e. □ 

Note that v E 'He may not lie in the domain of A''. If it did, one would have {A''u, v)e - 
{u,h''v)^, - {u,K''f)rw, where / - P^mV- The theorem could then be written S'q - 
U7=o ^''{f^in). However, this turns out to have contradictory implications. 

We now provide two results enabling one to recognize certain elements of S'q. 

Lemma 4.7. A linear functional f : Sq ^ C is an element of S'q if and only if there exists 
/? E N and Fq,F\, . . .F p e 'He such that 

p 

fiu) = 2<n, A'u)e, Vm e 'He- (4.7) 

Proof. By definition, f e S'q ifl 3p,C < oo for which |/(m)| < C||m||p for every u e Sq. 
Therefore, the linear functional 

ct> : dom(A*) ^ C by ^(u. Am, A^m, . . . A''m) = f(u) 

is continuous and Riesz's Lemma gives F - (/^jt)[^^o ^ ©Lo^^ with 
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P 

f{u) = (F, (u. Am, . . . A''m))0,^^ = Yj^P'k^ ° 
Corollary 4.8. If A .'He-^'He is bounded, then S'^ - 'He. 

Proof. We always have the inclusion "Kg ^ >Sq by taking p = 0. If A is bounded, then the 
adjoint A* is also bounded, and (4.7) gives 

f(u)^ly{A')'Ft,u\ , ^ueSa. (4.8) 

Since So is dense in "Kg by Lemma 4.2, we have / = Yjk=o(^*^''^k £ 'Ke- □ 

Remark 4.9. In view of Lemma 3.9, Corollary 4.8 shows that S'^ is a proper extension of 
"TYfi on any infinite network. 

In the case when the Gel'fand triple is constructed from the domain of A^^, as in Defi- 
nition 3.12, then one can extend A to distributions. 

Definition 4.10. Extend A to S'^ by defining 

mv,) := {6,,^hv, (4.9) 
so that A^(Vj:) - Y^y^ J Cvv(^(v^) - ^(vy)) follows readily from Lemma 3.16. 

S' ~ 

Now extend A to A defined on v_, e L^{^, P^) by A(vjc)(^) := A\\{^), so that 

(4.10) 

y~x 

Since v ^ i-> is an isometry, it is no great surprise that 

<v,,Av,,)i2= r vA^)U^)dFQ(0^{v„Av,)E. (4.11) 

5. The Wiener embedding and the space S'^ 

We have now obtained a Gel'fand triple So £ 'J^s £ (from either Lemma 3.10 
or Theorem 4.3), and we are ready to apply the Minlos Theorem to a particularly lovely 
positive definite function on So, in order that we may obtain a particularly nice measure 
on S'q. This allows us to realize bdG as a subset of S'^. Recall that So contains the energy 
kernel; see Lemma 3.11 or Lemma 3.15. 

5.1. The Wiener embedding. In [JP09c, §5], we constructed 'He from the resistance met- 
ric by making use of negative definite functions. We now apply this to a famous result of 
Schoenberg which may be found in [BCR84, SW49]. 

Theorem 5.1 (Schoenberg). Let X be a set and let Q : X x X —> R be a function. Then the 
following are equivalent. 

(1) Q is negative definite. 

(2) Vf e R^, the function pi(x,y) :— e"'2('^'>'' is positive definite on X xX. 

(3) There exists a Hilbert space "TY and a function f : X ^ such that Qix,y) — 
\\fix)-f(y)\t. 
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In the proof of the following theorem, we apply Schoenberg's Theorem with t - j to 
the resistance metric in the form 

R^(x,y)^\\v,-Vy\\l, (5.1) 

which appears in [JP09c, Thm. 2.13]. The proof of Theorem 5.2 also uses the notation 
Ef(/):= J^,J{^)dm). 

Theorem 5.2 (Wiener embedding). The Wiener transform 'W : "TYg — > L?(S'q, P) by 

'W-v^v, v(f):=<v,^>^, (5.2) 

is an isometry. The extended reproducing kernel [Vx]xe(? is a system of Gaussian random 
variables which gives the resistance distance by 

R^{x,y) = ^^{{v.-vy)\ (5.3) 
Moreover, for any u, v e "Hg, the energy inner product extends directly as 

{u,v)e, = (Sv) = ^ avdf. (5.4) 

Proof. Since R''{x,y) is negative semidefinite (see [JP09c, Thm. 5.4]), we may apply 
Schoenberg's theorem and deduce that exp(-^||M - v||^) is a positive definite function on 
'HgX'Hs. Consequently, an application of the Minlos correspondence to the Gel'fand triple 
established in Lemma 4.2 yields a Gaussian probability measure P on S'^^. 
Moreover, (2.19) gives 

E^(ei<".f>'H') ^ e-{Ml^ (5.5) 

whence one computes 



1 + i(M, - ^(m, + ■ ■ ■ j d¥(^) = 1 - i<M, M>£ + ■ • • . (5.6) 



Now it follows that E(m^) - E^{{u,^)^^,) - ||m||| for every u e Sa, by comparing the terms 
of (5.6) which are quadratic in u. Therefore, 'W : "Hg S'q is an isometry, and (5.6) gives 

(|vv - vvh = ((V, - v,„ ^f) = ||y, - y,|||, (5.7) 

whence (5.3) follows from (5.1). Note that by comparing the Unear terms, (5.6) implies 
Ef (1) = 1, so that P is a probability measure, and E^ ((m, ^)) - and E^((m, ^)^) - \\u\\^y, so 
that P is actually Gaussian. 

Finally, use polarization to compute 



<M,v)fi= ^ (||m + v|l| - 11m - v|||) 

= i (Ef (|fi + vf) - Ef (|S - vf)) by (5.7) 



I f |M + rf (f)-|M-rf (^)£/P(^) 
4 Js' 
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1 This estabUshes (5.4) and completes the proof. □ 

2 It is important to note that since the Wiener transform : So ^ S'q is an isometry, 

3 the conclusion of Minlos' theorem is stronger than usual: the isometry allows the energy 

4 inner product to be extended isometrically to a pairing on "He x S'q instead of just SgxS'q. 

5 Remark 5.3. With the embedding 'Kg — > L~{S'f^, P), we obtain a maximal abelian algebra 

6 of Hermitian multiplication operators L°°(>Sq) acting on L^{S'g,F). For a sharp contrast, 

7 note that the multiplication operators on f-ls are trivial, by [JP09b, Lem. 2.3]. This result 

8 states that if i/j : G" ^ R and denotes the multiplication operator defined by (M^m)(x) - 

9 (p{x)u{x), then is Hermitian if and only if = kl, for some ^ 6 R. 

10 Remark 5.4. The reader will note that we have taken pains to keep everything R-valued in 

11 this section (especially the elements of Sa and S'fj), primarily to ensure the convergence of 

12 J^, ei<">f>'v*' dF(^) in (5.5). However, now that we have established the fundamental identity 

13 {ii, v)e = Mvc/P in (5.4) and extended the pairing (■, to "He x S'q, we are at liberty 

14 to complexify our results via the standard decomposition into real and complex parts: u - 

15 Ml + iu2 with Ui R-valued elements of 'He, etc. 

16 Remark 5.5. The polynomials are dense in L^{S'q, P). More precisely, if ip{t\,t2, . . . , fj:) is 

17 an ordinary polynomial in k variables, then 

<P(^) -.^ ip[{Uu^)'W,{ul,0'W, ■ ■ ■ {Uk,^)'w) (5.8) 

18 is a polynomial on S'q and 

roly„ := {if {uiia U2(a ■ ■ ■ Uk(^)) , deg(^) < n, ■. uj e <Hfi, ^ 6 S'^] (5.9) 

19 is the collection of polynomials of degree at most n, and {Poly ^^Y^^^ is an increasing family 

20 whose union is all of S'q. One can see that the monomials {u,^)^, are in L^(>Sq,P) as 

21 follows: compare like powers of u from either side of (5.6) to see thatEf (<M,^)2«+i) = 

22 and 

e4<m,^)2«) = ^ \{u,0-^,^-" dm = , (5.10) 

23 and then apply the Schwarz inequality. 

24 To see why the polynomials {'Po/y„)J^Q should be dense in I}{S'q,¥) observe that the 

25 sequence {^'poh,, I^JIq orthogonal projections increases to the identity, and therefore, 

26 {Ppo/v„"} forms a martingale, for any u € 'Kg (i.e., for any u € L}{S'q, P)). 

27 Denote the "multiple Wiener integral of degree «" by 

H„ :- (cZspan{<M, ■)"^ : u € "He}) Q [{u, ■>^, -k <n,ue "He), 

28 for each n > 1, and Hq := CI for a vector 1 with ||1||2 = 1. Then we have an orthogonal 

29 decomposition of the Hilbert space 



03 

L\S'q,¥) ^ ^ H„. 
«=o 



(5.11) 
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See [HidSO, Thm. 4.1] for a more extensive discussion. 

A physicist would call (5. 1 1) the Fock space representation of L"(>Sq, P) with "vacuum 
vector" 1; note that H„ has a natural (symmetric) tensor product structure. Familiarity 
with these ideas is not necessary for the sequel, but the decomposition (5.1 1) is helpful for 
understanding two key things: 

(i) The Wiener isometry : He — * L-(S'q,¥) identifies TYe with the subspace Hi of 
L^(S'q,¥), in particular, L-{Sq,¥) is not isomorphic to TYg. In fact, it is the second 
quantization of He- 

(ii) The constant function 1 is an element of L^(S'q, P) but does not correspond to any 
element of '7/g. In particular, constant functions in "He are equivalent to 0, but this is 
not true mL}{S'a,V). 

It is somewhat ironic that we began this story by removing the constants (via the intro- 
duction of 6), only to reintroduce them with a certain amount of eff'ort, much later Item 
(ii) explains why it is not nonsense to write things like P(.Sg) - f , \dV — \, and will be 

helpful when discussing boundary elements in §6.1. 

Corollary 5.6. For e.v(^) := one has E^(e.v) = ^-^(o,.^) hence 

Ef(^ev)= r 7Ji)e,{0dV^e-^-^'^''y\ (5.12) 

Proof. Substitute u - Vj or m = Vj: - v,. in (5.5) and apply (5.1). □ 

Remark 5.7. Free resistance is interpreted as the reciprocal of an integral over a path space 
in [JP09c, Rem. 3.15]; Corollary 5.6 provides a variation on this theme: 

iJ^(x,y) = -21ogE^(F^e,) = -21og r JMeAOd'P. (5.13) 

Observe that Theorem 5.2 was carried out for the free resistance, but all the arguments 
go through equally well for the wired resistance; note that R^ is similarly negative semi- 
definite by Theorem 5.1 and [JP09c, Cor 5.5]. Thus, there is a corresponding Wiener 
transform : Tin L^iS'^, P) defined by 

: V /, / = P^,„ V and /(^) = (/, 0^y. (5.14) 

Again, {/i ) rec" is a system of Gaussian random variables which gives the wired resistance 
distance by R^{x,y) = E^((/, - fyf). 

Remark 5.8. For u 6 l-fann and ^ e S'^, let us abuse notation and write u for u. That is, 
m(^) := m(^) - {u,^}'ui. Unnecessary tildes obscure the presentation and the similarities to 
the Poisson kernel in §6. 

Remark 5.9. Theorem 5.2 showed that {s^} forms a system of Gaussian random variables. 
Since the Wiener transform is an isometry, 

Eie^) = and E(g,ev) = <^.t,v- (5.15) 

Since independence of Gaussian random variables is determined by the first two moments, 
it follows that {sj) forms a system of i.i.d. Gaussian random variables with mean and 
variance 1. This is noteworthy because while independence implies orthogonality, the 
converse does not hold without the additional hypothesis that the distributions be Gaussian. 
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6. The resistance boundary of a transient network 

With the tools developed in §3 and §5, we now construct the resistance boundary bdG 
as a set of equivalence classes of infinite paths. Recall that we began with a comparison 
of the Poisson boundary representation for bounded harmonic functions with the boundary 
sum representation recalled in (1 .2): 

u(x) — I u(y)k(x,dy) <-> u(x) — y + u(o). 



In this section, we replace the sum with an integral and complete the parallel. 

Corollary 6.1 (Boundary integral representation for harmonic functions). 

For any u € "Harm and with — P-HmmVx, 

«(x)= Lu{^)hA^)d¥^{^) + u{o). (6.1) 
Proof. Starting with (2.6), compute 



u{x) — u{o) — {hx, u)e — {u, hx}s — 




(6.2) 



where the last equality comes by substituting v = hx in (5.4). It is shown in [JP09b, 
Lem. 2.24] that JTy - hx. □ 

Remark 6.2 (A Hilbert space interpretation of bd G). In view of Corollary 6. 1, we are now 
able to "catch" the boundary between Sq and S'^ by using A and its adjoint. The boundary 

S' 

of G may be thought of as (a possibly proper subset of) Corollary 6.1 suggests that 
^{x,d^) :- hx{^)dF^ is the discrete analogue in 'He of the Poisson kernel k{x,dy), and 
comparison of (1.2) with (6. 1) gives a way of understanding a boundary integral as a limit 
of Riemann sums: 

I uhxdF^ = lim ^ u{x)'^(x). (6.3) 

(We continue to omit the tildes as in Remark 5.8.) By a theorem of Nelson, is fully 
supported on those functions which are Holder-continuous with exponent a = j, which 
we denote by Lip(i) c S'^; see [Nel64]. Recall from [JP09c, Cor. 2. 16] that ffe c Lip{\). 
See [ Arv76a, Arv76b, Min63, Nel69] . 

6.1. The boundary as equivalence classes of paths. We are finally able to give a concrete 
representation of elements of the boundary. We continue to use the measure P'^ from 
Theorem 5.2. Recall the Fock space representation of L^{S'q, P) discussed in Remark 5.5: 



where "Hg" :- CI for a unit "vacuum" vector 1 corresponding to the constant function, 
and "TYg" denotes the n-fold symmetric tensor product of "Kg with itself. Observe that 1 is 

S' 

orthogonal to Tin and 'Harm, but is not the zero element of L^(^, P'^). 
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Lemma 6.3. For all v e "Harm, ^s'^ v dF^ - 0. 



Tin 



Proof. The integral Jij, vd¥^ = Js'^ Ivt/P^ is the inner product of two elements in L^(^,P^ 

Tin 

which he in different (orthogonal) subspaces; see (5. 11). □ 

Alternatively, Lemma 6.3 holds because the expectation of every odd-power monomial 
vanishes by (5.6); see also (5.10) and the surrounding discussion of Remark 5.5. 
Recall that we abuse notation and write /i ^ = (/ii, = /i v for elements of S'q. 

Definition 6.4. Denote the measure appearing in Corollary 6. 1 by 

:= (1 + /i,.)c/P^- (6.5) 
The function 1 does not show up in (6. 1) because it is orthogonal to Harm. 

J^,^ u(l + h,) d¥^ = L u £/p2 + <M, = (m, /!^)£, for u e 'Harm, 

Till Till 

where we used Lemma 6.3. Nonetheless, its presence is necessary. 

Till Tin Tin Tin 

again by Lemma 6.3. 

Remark 6.5. We have shown that as a linear functional, ;U.i[l] = 1. It follows by standard 
functional analysis that fiy > P'^-a.e. on Thus, //^ is absolutely continuous with 
respect to P^ (yU j- «; P^) with Radon-Nikodym derivative ^ = 1 + h^-. 



Definition 6.6. Recall that a path in G is an infinite sequence of successively adjacent 
vertices. We say that a path oj = (xq, xi, X2, . . . ) is a path to infinity, and write w —> oo, iff 
a> eventually leaves any finite set F c G", i.e., 

3A^ such that n> N => x„ i F. (6.6) 

If = (xo,xi,X2, . . .) and W2 - {yo,y\,y2, ■■ ■) are two paths to infinity, define an 
equivalence relation by 

(Di ^ ll>2 lim(/!(x„) - h{y„)) = 0, for every h e M. (6.7) 

«— }oo 

In particular, all paths to infinity are equivalent when Harm - 0. 

If j6 = [<y] is any such equivalence class, pick any representative ai - {xq, xi,X2, ■ ■ .) 
and consider the associated sequence of measures {//.»„ ). These probability measures lie in 
the unit ball in the weak-* topology, so Alaoglu's theorem gives a weak-* limit 

vg := lim jdx,,. (6.8) 
For any h 6 Harm, this measure satisfies 

hd/jjc,, > hdvp. (6.9) 
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1 Thus, we define bd G to be the collection of all such yS, and extend harmonic functions to 

2 bdGvia 

h(/3):^ f hdvp. (6.10) 

JbdG 

3 Definition 6.7. For u e TYg, denote ||m||oo su^^^qo \ii{x) - u{o)\, and say u is bounded iff 

4 ||m||oo < OO. 

5 Lemma 6.8. Ifv e "Hg /i bounded, then PrmV is also bounded. 

6 Proof. Choose a representative for v with < v < /(T. Then by Corollary 6. 1 and (6.5), 

P-H..,Ax)^ r, v(f)/i,(^)Jpe(^) + M(o)= L^^i^)d^iA^) + u{o). 

7 Since yu.v is a probability measure (cf. Remark 6.5), we have P-Hm-mV ^ 0, and hence the 

8 finitely supported component Pjri„v - v — P.HcrmV is also bounded. □ 

9 Lemma 6.9. Every v e Mis bounded. In particular, \\Vx\\oa < R^(o, x). 

10 Proof. According to Definition 2.15, it suffices to check that Vv, w\. and are bounded 

11 for each x. Furthermore, ||w^v||oo = ll^'rmM'vll'» - ll*^tll'» Lemma 6.8, and w\ = v^ + Wg by 

12 definition, so it suffices to check and Wg. By [Soa94, Lem. 3.70], w,, has a representative 

13 which is bounded, taking only values between and Wo{o) > 0. It remains only to check 

14 Vj. The following approach is taken from the "proof" of [JP09c, Conj. 3.18]. 

15 Fixx,y e G'^ and an exhaustion {G/.}, and suppose without loss of generality that o,x,y e 

16 Gi. Also, let us consider the representative of specified by Vx(o) = 0. On a finite 

17 network, it is well-known that 

Vx^R{o,x)Ux, (6.11) 

IB where Ux(y) is the probability that a random walker (RW) started at y reaches x before o, 

19 that is, Uxiy) :- ^y[Tx < To], where Tx denotes the hitting time of x. This idea is discussed 

20 in [DS84,LPW08,LP09]. 

21 Therefore, one can write (6.1 1) on Gk as v^*' = Rc^{o, x)uf\ In other words, v^P is the 

22 unique solution to Av - 6x - do on the finite subnetwork G^. Consequently, for every k 

23 we have vf\y) < Rqf{o,x) for all y 6 Gk- Since R^{x,y) - lim,t^oo ^y [JP09c, 

24 Def. 2.9], we have ||Vf||oo < R^{o, x) for every x e G°. □ 

25 Theorem 6.10. Let /3 ehdG and let to — (xq, xi,X2, . . .) is any representative of p. Then 

26 yS e bd G defines a continuous linear functional on So via 

/3(v) lim I vd^l,, v e So- (6.12) 

27 In fact, the action of p is equivalently given by 

Piv) = lim P^HanMXn) - PlianMo), V 6 Sa- (6.13) 
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Proof. To see that (6.12) and (6.13) are equivalent, compute 



Vhr 



{Vjlx„)e. - P'Ha„nV(Xn) - P'Ha,mV{0), 



\ y(l + /J.vjc/P^ = fcWf^ + r 

because 1 is orthogonal to TYe in L-{Sq, F^); see (5.11). 

Now, to see that (6.12) or (6.13) defines a bounded linear functional, we only need to 
check that sup„£^g{yS(v) i ||v||£ = 1) is bounded, but this is the content of Lemma 6.9. Note 
that the equivalence relation (6.7) ensures that the limit is independent of the choice of 
representative co. a 

Remark 6.11. In light of (6.13), one can think of vg in (6.8) as a Dirac mass. Thus,/? e bdG 
is a boundary point, and integrating a function / against corresponds to evaluation of / 
at that boundary point. 

7. Examples 

Our presentation of bd G may appear somewhat abstract in the general case. However, 
we now illustrate the concept with a simple and entirely explicit example where the rep- 
resentation by equivalence classes given at the end of §6 takes on an especially concrete 
and visual form. Moreover, the computations can be completed without the direct con- 
struction of Sg, Sq, or any discussion of Lr{S'(;, P); we can obtain the boundary simply by 
constructing certain functions on the network. We feel this is an especially nice feature of 
our approach. 

Example 7.1 (One-sided infinite ladder network). Consider two copies of the nearest- 
neighbour graph on the nonnegative integers Z^, one with vertices labelled by {x„), and the 
other with vertices labelled by {y„}. Fix two positive numbers a > 1 > j6 > 0. In addition 
to the edges Ct___v, j - a" and c,, , = a", we also add "rungs" to the ladder by defining 



Xq 



■ X\ 



X2 ■ 



■ Xi 



(7.1) 



p 






a 







P' 



yo y\ yi jn 

This network was suggested to us by Agelos Georgakopoulos as an example of a one- 
ended network with nontrivial 'Harm. The function u constructed below is the first example 
of an explicitly computed nonconstant harmonic function of finite energy on a graph with 
one end (existence of such a phenomenon was proved in [CW92]). Numerical experiments 
indicate that this function is also bounded (and even that the sequences {u{x„)}'^^f^ and 
{m(3';i)))JLo actually converge very quickly), but we have not yet been able to prove this. 
Numerical evidence also suggests that A is not essentially self-adjoint on this network, but 
we have not yet proved this, either. 

This graph clearly has one end. We wiU show that such a network has nontrivial resis- 
tance boundary if and only if a > 1 and in this case, the boundary consists of one point for 
/3 - 1, and two points for /3 such that (1 + ^)^ < a/fS'^. It will be made clear that the paths 
Wi - ix\,X2, X3, . . . ) and ojy - (yi,3'2,3'3, ■ ■ ■ ) are equivalent in the sense of Definition 6.6 
if and only if yS = 1 . 

For presenting the construction of u, choose /? < 1 satisfying 4ff < a (at the end of the 
construction, we explain how to adapt the proof for the less restrictive condition (1 + ^)^ < 
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a Iff). We now construct a nonconstant u 6 'Harm with m(jco) = and uiya) - -1. If we 
consider the flow induced by m, the amount of current flowing through one edge determines 
u completely (up to a constant). Once it is clear that there are two boundary points in this 
case, it is clear that specifying the value of u at one (and grounding the other) determines 
M completely. 

Due to the symmetry of the graph, we may abuse notation and write n for x„ or y„, 
and « for the vertex "across the rung" from n. For a function u on the ladder, denote the 
horizontal increments and the vertical increments by 

du{n) :— u(n + 1) - u{n) and cruin) :— uin) - u{n), 
respectively. Thus, for « > 1, we can express the equation Am(«) = by 



Am(«) = a"5u{n - 1) - a"^^6u{n) + /3"o-u(n) = 0, 
which is equivalent to 

1 P" 
6u(n) — —6u(n - 1) H -cru{n). 

a a"+' 

Since symmetry allows one to assume that u(n) = 1 - u(n), we may replace cru{n) by 
2u{n) + 1 and obtain that any u satisfying 



„(„ + 1) . „(„) + ilitti) + 1 «(„) + i (^)" (7.2) 

is harmonic. It remains to see that u has finite energy. 

Our estimate for &{u) < oo requires the assumption that a > 4/3^, but numerical compu- 
tations indicate that u defined by (7.2) will be both bounded and of finite energy, for any 
/3 < I < a. First, note that m(1) - ^ and so an immediate induction using (7.2) shows that 
6u{n) - u{n + 1) - u{n) > for all n > 1, and so u is strictly increasing. Since /3 < I < a, 
we may choose so that 



n>N 




Then n > N implies 

M(n + 1) < 2m(«) + -, (7.3) 
a 

by using (7.2) and the fact that u(n) is increasing and ^ < 1- Now use (7.2) to write 
6u(n) = ^(6u)(n - 1) + m(«) + i) (f )" 

n-I 

= ^(5m)(0) + 2?(^(«-^> + ;j)(S"' 

where the second line comes by iterating the first, and the third by algebraic simplification. 
Applying the estimate (7.3) gives 
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2 ^/M(fc) < 2^ Y^p'uik - 1) + 1 ^/?* = 2^ Y^p'uik - 1) + 2f 

k=l k=\ k=\ k=2 

and iterating gives 



i-p 



^ ^ p{\-n ^ (2/3)" ^ y ^jP" - P" 



-P 



(7.4) 



Now the energy &{u) - 2,^o {6u{n))^ can be estimated by using (7.4) as follows: 



o(m) < > r 1 H ^ 

/7=0 ^ ' 



Pi I-P") (2p)" 2p + 2/3"+ ' - 2"+2y3«+ 1 - 22j6«+2 + (2;8)"+2 \^ 



a(l-p)(2p-l) 



and the condition a > 4/?^ ensures convergence. 

Note that this computations above can be slightly refined: instead of a > 4/?^, one need 
only assume that a > (I + ^)^P^- Then, fix e > for which ajp" > (1 + ^)^ + e and choose 

so that n > N implies (PI a)" < l + ^+e(l+2Q; + as). Then the calculations can be 
repeated, with most occurrences of 2 replaced by 1 + ^ + e. 

Remark 7.2. [Comparison of Example 7. 1 to the 1 -dimensional integer lattice] In [JP09b, 
Ex. 6.3], we showed that the "nonnegative geometric integers" network 



supports a monopole but not a harmonic function of finite energy, for ff > 1 . These con- 
ductances correspond to the biased random walk where, at each vertex, the walker has 
transition probabilities 



p{n, m) - 




, m — n — I, 
, m — n + I. 



In particular, this is a spatially homogeneous distribution. In contrast, the random walk 
corresponding to Example 7.1 has transition probabiUties 



1 

l+ff+(f)"' 


}n 


— n 


- 1, 


l+ff+(f)"' 


m 


— n 


+ 1, 


(J3/a)" 


m 


— h 











Thus, Example 7.1 is asymptotic to the nonnegative geometric integers. 

One can even think of Example 7. 1 as describing the scattering theory of the geometric 
half-integer model, in the sense of [LP89]. In this theory, a wave (described by a function) 
travels towards an obstacle. After the wave collides with the obstacle, the original function 
is transformed (via the "scattering operator") and the resulting wave travels away from the 
obstacle. The scattering is typically localized in some sense, corresponding to the location 
of the collision. 
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1 To see the analogy with the present scenario, consider the current flow defined by the 

2 harmonic function u constructed in Example 7.1, i.e., induced by Ohm's law: I(x,y) - 

3 c„(m(x) - u(y)). With drv.Jx) := ^ 2 \I{x,z)\, this current defines a Markov pro- 

" ^ iz : /(A-,z)>oi 

4 cess with transition probabilities 

diV|,|(jc) 

5 and P(x,y) = otherwise; see [JP09d, JP09a]. This describes a random walk where a 

6 walker started on the bottom edge of the ladder will tend to step leftwards, but with a 

7 geometrically increasing probability of stepping to the upper edge, and then walking right- 

8 wards off towards infinity. The walker corresponds to the wave, which is scattered as it 

9 approaches the geometrically localized obstacle at the origin. 
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